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Abstract 

C . 

| The Vlasov type quantum kinetic equation for the deconfined quarks in the 

00 ■ strong quasi-classical gluon field is derived in the covariant single-time formalism. 

The equations system for the Wigner function components is obtained as a result 

. of spinor and color decomposition. The field-free and vacuum solutions of the ki- 

netic equation are found, and the conservation laws are derived. The flux-tube 

\ configuration of gluon field is discussed in detail. 

^ : 1 Introduction 

o ; 

The creation and evolution of the quark-gluon plasma (QGP) by ultra-relativistic heavy 
ion collisions is a very complex process occurring on the time scale of lfm/c and at 
energy density above of lGev/ fm 3 . Parton gas probably comes to equilibrium state 
• in the process of hadronization. Details of this transition are of great importance for 

>• ■ interpretation of experimental data. If the QGP is really formed in experiments at the 

^ . RHIC, its thermalization time must be so small, that it presents a serious problem for 

theoretical explanation pQ. Several theoretical tools are applied to describe the various 
physical phenomena accompanying the collision. The kinetic equation is the basic means 
for investigation of the nonequilibrium evolution PJ-P-7 . The flux tube model JB] and 
the Schwinger mechanism of pair production are often used for the research of an early 
stage of the QGP formation. Vacuum pair creation is on essentially non-perturbative 
effect, which description requires the exact solution of the field equations. The resulting 
source term in the kinetic equation (KE) describes the production rate and momentum 
spectra of the created particles [T9] . 

The QED example shows that the evolution of a plasma created from vacuum is rather 
dependent on the structure of the source term [213- [22] • The semi-phenomenological 
source term based on the classical Schwinger formula (e.g., [231211) can result in signif- 
icant inaccuracy for the fast varying electrical field. In particular, such source can not 
reproduce correctly the relation between the production rates of two components with 
different masses and statistics EH E3D • The main reason is that correct source term has a 
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non-markovian time dependence ^3 120] , as against the Schwinger formula. The most in- 
teresting effects arising from these feature are the suppression of boson creation with zero 
kinetic momentum and the suppression of statistical factor influence. The Schwinger- like 
source contains the statistical factor as multiplier and produce the suppression of fermion 
creation (Pauli blocking) and enhancement of boson creation. The correct source term 
[2*Uj contains the statistical factor as integrand on time (non-markovian property) that 
facilitates the fermion creation and reduces the boson creation by large particle density. 
The joint action of these factors may cause the effect of "statistic inversion" at the short 
times scale when the fermions production rate is more greater than bosons one. This 
effect can be so strong that heavy fermions are created more actively than light bosons. 

The other noteworthy feature of a proper source term is that a momentum distribution 
of the produced fermion pairs is close to the quasi-equilibrium one both in the transverse 
and in longitudinal directions. The influence of such source term can facilitate more 
rapid formation of a quasi-equilibrium state of QGP along with the small mean free 
path compared to the Compton length. The correct source term allows to separate the 
vacuum polarization effects which can exceed the contribution of real particles in the 
thermodynamical variables |2*7j . 

The kinetic equation with account of vacuum creation effect can be obtained in the 
Wigner function approach, which is widely used for the description of relativistic quantum 
systems during the past few decades. It is more convenient to use the single-time Wigner 
function variety [2B1 EHj to solve the initial value problem for the collisions at RHIC. A 
covariant version of this approach has been recently developed in jHO] for the QED case 
within the framework of the proper time method on space-like hyperplanes. We use this 
approach in the present work for the derivation of the kinetic equation for the covariant 
single-time Wigner function (WF) in the quark sector of QCD at presence of a strong 
quasi-classical background gluon field. 

The paper is organized as follows. We introduce the equations of motion for quarks and 
quasi-classical gluon fields in covariant proper time approach in Sect. 2. The single-time 
Wigner function is introduced here as the basic element of kinetic theory. The physical 
observables are bilinear compositions of quark field operators, and can be calculated via 
WF. The kinetic equation for the Wigner function is derived in a matrix form in the 
spinor and color spaces in Sect. 01 The spinor decomposition of this KE is also fulfilled 
for the more convenient analysis. This representation of KE is useful for construction of 
conservation laws, Sect. El The vacuum and field- free solutions are obtained in Sect. |H 
The vacuum WF play a role of initial value for a solution of the Cauchy problem. The 
special sample of the kinetic equation is investigated in Sect. El which is inspired by the 
flux-tube ^H] model of vacuum quark production under conditions of ultra-relativistic 
heavy ion collisions. Comparison with QED case is made in Sect. Finally, Sect. El 
summarizes some results of the work. We use the system of units with h — c — 1 and the 
signature of metric tensor (1,-1,-1,-1). 
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2 Basic Equations 

We start with the QCD Lagrange density (for one quark flavor, g > 0) 

Cqcd = - igA^ - mi>i> - \f^\ (1) 

with the quasi-classical gluon field 

A^ = A«t a , F lu , = d li A v -d v A li -ig[A li ,A v ], a=l,2...N, (2) 

where t a = X a /2 are generators of SU(N) gauge group in fundamental representation. In 
particular, A a are the Gell-Mann matrices for SU(3) and the Pauli ones for the SU(2) 
group. 

The covariant proper time derivative d T is defined by means of the space-time parametriza- 
tion by a family of space-like hyperplanes a(n, r) 

(7(n,r)±n M , n% = r, n 2 = I. (3) 

Hence there is a covariant decomposition 

x n _ rn M _|_ x ^ _ n ^ _|_ Q±^ 

d T = n^ d^ = A, u d\ bT = <r-n»n v , (4) 

where is a unit time-like vector, A^ denotes the transverse projector, and g^ v is the 
metric tensor. Analogous decomposition is produced for any vector / M 

/" = /"n" + /£, /»=/%, K = Kr, (5) 

and for any anti-symmetric tensor c Miy 

= fn* - cPr? + d? , c» = n u c^, cf = A£A£c Q/? , (6) 

so = c M . In these terms, the basic equations of motion for quark field operators are 

<9 T ^ = -t"7±(^ - igAp)tl> - im^ip + igA"ip, 

d T ip = -iify +igA$iW +1771^7" -ig^AK (7) 

The mean gluon fields obey the equations 



9^ = (8) 

where -E^ = n^F^ and represent the color electrical and color magnetic fields, re- 
spectively, and J u is a color current 

J V = JX, J2=<hvt a i>>a, (9) 
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the symbol < ... > CT designates the average with statistical operator of the system on a 
hyperplane in the Heisenberg picture. 

The basic element of statistical description is the covariant Wigner function (WF) on 
the space-like hyperplane er(r) |2H1 13U]: 

x U™'(x,x + y/2)pf(x,y,T)U b ' b (x-y/2,x), (10) 

here the upper and bottom indices are the color and spinor ones, respectively, p is the 
one-particle density matrix of quarks 
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Pik(x,y,r) = -- < 



^(x + y/2),iP b k (x-y/2) > CT , (11) 



and U is the unitary link operator 

U(xi,x 2 ) = exp I ig I dz^A^z) I , (12) 



provides the gauge invariance of the WF 29J. The integral is taking along straight line 
between points x\ and x 2 connected with space-like interval and 

U + (xi,x 2 ) = U~ 1 (x 1 ,x 2 ) = U(x 2 ,X!), 

U(xi,x 2 ) U (x 2 , X\) = U{x\,Xx) = 1. (13) 

The dynamical variables can be expressed in terms of the Wigner function. For ex- 
ample, the color and electromagnetic current densities are 

J*(x) = Tr < lu t a W>, (14) 
j v {x) = Tt< lv W>, (15) 

where trace Tr is carried out on spinor and color indices. The brackets < ... > denote the 
covariant momentum average on the hyperplane p ■ n = jJUJ, e.g. 

<W >= J -^^5(p-n)W(x ± ,p ± ;T = x-n). (16) 

The integral in (jlfjjl provides the locality of corresponding observables, due to that function 

8r(?) = J ^e-*-i(p.n) (17) 

plays the role of the three-dimensional delta function on a hyperplane a(n, r). 

The color WF have a rather complex matrix structure, therefore it is convenient to use 
the corresponding decompositions in the spinor and color space. The spinor decomposition 
in a complete basis of the Clifford algebra is 

W = al + + c^ v + rf M 7 M 7 5 + *e 7 5 , (18) 
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where the coefficient functions a (scalar), (vector), c^ u (antisymmetric tensor), d (axial 
vector), e (pseudo-scalar) are the hermitian color matrices 

a = -trW, & M = -tv(YW), cT = ^tr(a^W), 

d" = --tr(y7 5 W), e = ~tr (-y 5 W), (19) 

the symbol "tr" denotes the trace on spinor indices only, and 7 5 = — ry°7 2 7 2 7 3 . 

The algebraic structure of the WF in a color space of N x N matrices is represented 
by the observable color singlet W a and the unobservable color multiplet W a , 

W = W s -i + W a -t a , a=l,2...(iV 2 -l), (20) 
W s = (1/N) tr c W, W a = 2tr c (t a W), 

where 1 is unit matrix and tr r is the trace on color indices. 



3 Kinetic Equation 



We calculate the proper time derivative of Eq. ljlUJl for a derivation of the KE for the 
Wigner function 



d T w£ 



d 4 y e ip - y 6{y ■ n) {U aa '{x, x + y/2) U b '\x - y/2, x)d T pf k b ' + 
+Pfk'dr [U aa ' (x, x + y/2) U b '\x - y/2, x)] } 



:2ii 



and substitute the time derivative of the field operators from the field equations (JJJ). The 
emerging terms with perpendicular derivative d x are transformed via the integration by 
parts. The derivative rules for link operator are follow from known formula 12H] 

SU (xi, x 2 ) = ig8xiA fl (xi)U(x 1 , x 2 ) — igU (xi, x 2 )A^(x 2 )5x 2 — 

i 

-ig J dsU{x u z{s))F^{z{s))U{z{s), x 2 ){x 1 - x 2 )»[5x 2 + s{5 Xl - 5x 2 )] u , (22) 
o 



where z(s) = x 2 + s(xi — x 2 ), e.g. 

dJx)U(x, x + y/2) = -ig 



AW(x,x + y/2)-AJx) + 



1/2 

+ J dsFW(x + sy)y» 



U{x,x + y/2), 



(23) 



where the notation for the "Schwinger string" 

AWz) = U(x,z)A ll (z)U(z,x) 



(24) 
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is introduced. Then we express the variable y as the momentum derivative 

y^exp (ip ■ y)5(y ■ n) = 6(y ■ n)d^(p) exp (ip ± ■ y ± ), 
and we obtain after some algebra the exact equation of motion for the WF 



(25) 



d T W + -d^(x)[S», W] + imfr 11 , W\ - ip^{S", W} + 

1/2 

+g&t(p) J ds {eP(t,x ± - isd ± (p))W + WEW(t,x ± + isd ± (j>)) + 
o 

+^?J (r, x ± - isd ± (p)) ([W, S»] - 2s{W, + 
+\(\W i S v ] + 28{W t S v }yS! (r,^+isd ± (p))} = 
= ig[Ai(x),W] + %[S'>,[Ai,W\ 



where 



gn = 7 ii 7 m 



AJJ7". 



(26) 



(27) 



This equation is resulted in a local form by means of the gradient expansion of a gluon 
field. We obtain in the first oder of this procedure the kinetic equation of the Vlasov type, 
which is correct for a rather slowly changing gluon field 

d T W+ l -d^x)[S^W} + 
+ 9 ~dl(p) (4{W, E,} + 2 {F, v , [W, S»]} - [F,*, W S"})) = 
= ip^S", W} - im[yKW] + ig ([A {l (x), W] + [A^ , [S», W}]) . (28) 

This equation despite of rather compact form is a very complicated matrix one in the 
direct production of 4 x 4 spinor and N x N color space. It is convenient to expand the 
WF in a some basis of this space for the separation of a different physical contributions. 

We perform at first the spinor decomposition (|T5j) in the KE. We obtain the equations 
set calculating the traces of Eq. ()28j) with the basic matrices of Clifford algebra 



(29) 



d T a + |<9» ({a, E a } + i[F^, c?]) = 4p% + ig[A\a], 

d T b" + n"<9X - <9^H + 
+§??(P) ({&", E a } + g^{b\F^} - n^, F£p} + % -[F^ ds\et*» 

= 2ptd^ + AmJt + ig ([A\W] + n"[A°, &£] - [A*, , (30) 
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+f d?(p) ({E a , cT} + {F^, (n^cf - n v df + g^c v - g^c»)} + 



2 

i 

+ 4 



+^L4", c""] + ig [A^, (n"c^ - n"c^ + ^c" - ^c")] , (31) 

d T d"-d^ + n^X + 

= 2j£b ft e af, ' i + 2mn"e + ^ d"] + r^L4?, c£] - K, d"]) , (32) 

9 r e + |9J(p) ({e, £ a } - |[F a/3 , c p ]e^ 

= -2p QC/35 e aps - 2md II + ig[A\ e], (33) 

where £ a/31 denote the convolution of the normal n M with the totally anti-symmetric unit 
tensor £ a ^ 6 \ i.e. 

e a ^ = n p e pa ^, e 0123 = +1. (34) 

It is more convenient for the analysis of concrete field configurations to rewrite that system 
concerning a projections on time-like and space-like directions 

d T a + §d» ({a, E a } + i[F^ <£]) = 4p% + ig[A\a], (35) 
d T bW + 9;(x)6i + §<9» ({6", - {fe^, F^}) = ig [A a , b a ], (36) 

drK + |9?(p) ({&?, + {b\F^} + l -[F^ dj]e^ 
= d^(xp + 2p^e a(3fM + Amc» + ig ([A\bf\ - K, fe 11 ]) , (37) 



9 T c^ + dt{x)c7 + ±8fXp) ({E a , c»} - {F^ } + a] 



= - mfe^ + i 5 ([A", C y + c^]) , (38) 
9 T cf + dr ± (x)c*-&t(x)J' ± + 

= epi + ig ([A", cf ] + K, c£] - [Af <]) , (39) 

+ 0?(*R + l&Hp) {{d\E a } - d{}) = 

= 2me + ig([A^S] + [A a ± ,d%) , (40) 
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flrdf + f 9?(p) ({<> + <*"} + ~[^, ttle*") = 

= + JZp^e *" + ig ([A\ d*\ - K, d"]) , (41) 

d T e + §<9» ({e, E a } - % -[F^ cj p )e^ = 

= -2pt4 s e^ s - 2md*+ig[A\e]. (42) 

The non-abelian character of gluon field is displayed in particular in the presence of the 
vector potential in right-hand side of these KEs that is necessary for the gauge invariance 
of the theory. 

It is the resulting set of KE for the description of deconfined quarks in a strong quasi- 
classical gluon field. 



4 Field-Free Limit and Vacuum Solution 

The equations system (|35j) - (j42j) can be solved exactly in the field-free limit = 0. We 
obtain assuming that all derivatives vanish: 

c » = 0> C U = °> d » = °> e = o> 

p^a(p) - mb^(p) = 0, 

p^e a ^ = 0. (43) 

The general solution of (J4HJ) is 

W( v )= v ^l (44) 
m 

with an arbitrary momentum dependence of a(p). The equilibrium WF [33] or the vacuum 
solution follow from the Eq. (pf4*|) as a particular cases. We perform the direct calculation 
of the function (fTUj) to select the vacuum case in the solutions class (02). The fermion 
operators on the hyperplane obey free anticommutation relations [3UJ EI] 

{^(x 2 ),V b (xi)} = 5 ab ^5 T {x x -x 2 ). (45) 

Then allows to write down the commutator (fTT]) as 

p ab (x,y;r) =<: <(z 2 )^ k (zi) :> T - 

^ / TTTT^ ^<5(£; 2 - m 2 )(m - jfe • 7 ) ^(A; 11 ) + 0(-fcH] , (46) 
2 7 (27r) d 

where the symbol : : indicates normal ordering. We obtain the vacuum WF by carrying 
out the averaging on vacuum state and substituting this representation in Eq. (fTU|) 



where u)(p±) = y m 2 — p\. This function is degenerated in the color space as a consequence 
of the primary supposition about the Lagrange density structure (JTJ). The WF (J4TJ) is the 
base for a solution of the Cauchy problem for the equations system (I35J) - (J4*2"|) . It can be 
proved that the vacuum solution, as well as the free-field ones, gives no contribution to 
the current densities (fT^) . (fTS|) . 

The kinetic equation has in the field-free case besides the vacuum solution (J4TJ) the 
another solutions of the type 

W a (p) =r] a W s (p) } a = l,2...N, W = {a, b^ . . . 0}, (48) 

where if is an arbitrary real numbers. 

5 Conservation Laws 

We calculate the divergence of electromagnetic current (|T3j) . using the spinor and color 
decompositions (JTHJ) and (|2T)j) 

d^(x)j„(x) = ANd^ix) < b^x) >= AN[d T < b\ > +d»< b^ s >]. (49) 

Performing the momentum averaging procedure (|16J) and taking the trace in Eq. (J36J) . we 
obtain the electromagnetic current conservation law 

d^(x)j^x) = 0. (50) 

We multiply the equation ()36|) on matrix t a and repeat the same procedure to derive the 
equation for the color current J® = (1/2) < b a ^ >. The result is 

d»(x) J*(x) + gf abc A b ^ c = 0. (51) 

The energy density e q of quark matter corresponding Eq.(jT|) is 

e q = n^n u T^ = i : ^{d T - igA*)ij) :, (52) 

where T£ v - energy-momentum tensor. Using the field equations (JTJ), we have 

e q =: i'f r y' L (d ± — igA ± )ip + rwi/ji/; :, (53) 

This variable can be written in terms of the WF as 

e g = Tr < (m - j L p x )W > +e vac , (54) 

where e vac = 2N < ui(p ± ) > is divergent vacuum contribution. We obtain using the spinor 
and color decompositions 

E M = 4N[m < a s > - < p^ s >] + E vac . (55) 

The linear combination of Eqs.()35p and (J3TJ) is used to calculate the right side of this 
equation. The result is 

d T E M + 4d*(x) < p^b\ > +gJ^ a ET = 0. (56) 
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6 Space-Homogeneous Color Field 



6.1 Flux-tube field configuration 

We consider here the "instant" frame of reference where n M = (1,0,0,0) and the field 
configuration typical for the flux tube model 

A; = (0, 0, 0, A a (t)), F % = -A a 3 (t) = E a 3 (t) = E a (t), (57) 

where the dot denotes the derivative with respect to time r = t and the Hamilton gauge 
was selected. We have E = (0, 0, E), H = in 3- vector representation 

F>» = (-E, H), <T = (-c 1; c 2 ), c" = (0, - Cl ). (58) 

We limit oneself to simple sample of SU(2) group (a = 1, 2, 3) below, where 

{t a ,t b } = V b i, [t a ,t b ] = if abc t c , ti c (t a t b ) = h ab , (59) 

the totally anti-symmetric structure constants f abc coincide here with the totally anti- 
symmetric unit 3-tensor e ljfc , f 123 = +1. As a result of color decomposition of the WF, 
the system (jHHJ)-(JI2J) is reduced to (d p = d/dp 3 ) 



d t a s + ^E a d p a a 


= 4pc*, 




d t a a + gE a d p a s 


= 4pc?, 


d t bl + 9 -E a d p % 


= o, 




d t b a +gE a d p b s 


-gf abc A b b c , 


d t b s + ^E a d p h a 


= 2p x 


d s + Amc\, 


d t b a + gE a d p b s 


= 2p x 


d a - Amc* + gf 


d t c{ + |£ a <9 p c? 


= -pa s 


+ mb s , 


dtcl + gE a 8 p c{ 


= -pa a 


+ mb a + gf ab % 


dtC s + g -E a d p c a 2 


= -pe s 




d t c a 2 + gE a d p c s 2 


= — pe a 


+ gf abc (A b x cl 


d t d s + 9 -E a d p d% 


= 2me s 




d t d a +gE a d p d s 


= 2me a 


+ gf abc A b c{, 


d t d s + ^E a d p d a 


= 2p x 


b s , 


d t d a + gE a d p d s 


= 2p x 


h a + g fabc A b d c^ 


d t e s + ^E a d p e a 


= 4pc* 


- 2md s , 


d t e a + ^E a d p e s 


= 4 P c« 


- 2md%. 



abc a b\c. 
A Uq, 



(60) 
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The important difference of these equations from their QED analogue is the structure of 
a force terms: the derivatives on time and on momentum act on different parts of WF 
(singlet and multiplet, respectively). This feature does not allow to reduce the problem 
to solving of ordinary differential equations even for the simple field configuration (|57|). 
The situation becomes complicated even more in case of SU(3) (a = 1, 2, . . . , 8) where 

{t a ,t b } = \5 ab l + d abc t c , (61) 

and d abc are the totally symmetric structure constants. We shall write out for an illustra- 
tion two first the equations (JoTij) only: 

d t a s + -E a d p a a = 4pc s 1; 
6 

d t a a + gE a d p a s + gd abc E b d p a c = 4pc?. (62) 



6.2 Constant Chromo-Magnetic Field 

Now we consider the gluon field configuration corresponding to a space-time homogeneous 
field 

A» = (0, A), = (0, H), A = const, H = const. (63) 

In this case the system (|3~5j) - (|32~j) is reduced to two independent groups of the equations 
(D = ?Hx d p ) 

{D,b} = 2ig[A,b], 
z[D, c 2 ] = -8p c 2 + 4md° . 

{Dx,c 2 } + -[D,a] = 2 (pa - mb + ig[Ax, c 2 ]) , 

{D,d°} + ^[Dx,b] = -4p x b-2ia[A,d°], (64) 

{D,d} = -4me + 2i#[A,d], 
z[D,ci] = -8pci, 
-i[Dx,d] = 8p x d+ 16mci — Aig[A, b°], 

{Dx, Cl } + ^[D,e] = 2(-pe + i< 7 [Ax,c 1 ]). (65) 

The solution for a corresponding abelian QED case is known [2Bj, but the derivation of 
a corresponding non-abelian analog is rather time-consuming work. We are limited here 
with a very simple case that have no analogy in QED. Assuming that A a = A, Va ("colour 
democracy"), we have H = and the system (fBl|) .([fiH j) is reduced to 

[A, b°) = 0, [A, b] = 0, pa = mb, (66) 

all other components are zero. This system allows the solutions with zero colour currents 
but with non-zero colour charges, for example 

b° ~ a, a = ail + a 2 A, b = ap/m, (67) 

where ai,a 2 are arbitrary constants. 
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6.3 Chiral limit 

The additional simplification is possible in the chiral limit m — > 0. In this case, the 
equations for the vector components of the Wigner function are separated from others. 
We suppose also the flux-tube symmetry in the electric field direction n = ~E/E, then 

b = bn + b ± p ± , 

d = d-(nxp ± ), (68) 

because of b is polar and d axial vectors. By the initial conditions of the type (|48jl . it 
follows from Eq. fjfiOj) that d° = and b° = (neutral system). Then we obtain 

d t b s + ^E a d p h a = 2pxd s 

d t b a + gE a d p b s = 2p x d a , 

d t d s + ^E a d p d a = 2p x b s , 

d t d a + gE a d p , d s = 2p x b a , 

fabc A b h c = Qj 

fabc A b d c = Q (QQj 

The two last algebraic equations play the role of constraints and have the particular 
solution ("color democracy") 

A a = A, b a = b, d a = d, a = 1,2, 3. (70) 

These conditions correspond to the special case of Abelian dominance approximation [HHU 
in relation to the gluon field. If we assume that the conditions (fTUj) are carried out for all 
others components WF also, the solution get the form (jlH|) 

W a (t,p)= V W s (t,p), a = 1,2,3, (71) 

where r\ is some parameter. We find two admissible values after substituting that in the 
system (JBUj) 

" = ± ^ (T2) 

The system (ffiTIj) is reduced to the three scalar equations at the account of the represen- 
tation flEEl 

d t \ + 7] gEd p fey = 2p\d, 
d t b ± + r]gEd p b ± = -2p {l d, 
d t d + V gEd p d = 2(pn6 x -6,|). (73) 

This equation set can be solved numerically by the characteristics method. The remaining 
part of Eq. (JoTij) is reduced to 

d t a + r]gEd p a = pci, 
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d t ci + rjgEdp Ci = -pa, 
d t c 2 + r]gEd p C2 = -pe, 
d t e + rjgEdp e = 4pc 2 . (74) 

We find two integrals of motion combining these equations in pairs 

D t (a 2 + Ac\) = 0, -> a 2 + 4c 2 = const, 

A(e 2 + 4c 2 ) = 0, -> e 2 + 4c 2 , = const, (75) 

here -D t = d t + i] gEd p . It follows for the field-free initial conditions that e = and Ci = 0. 
The self-consistent evolution of the mean gluon field obeys the Yang-Mills equation 

The equations (|73p. (|74|) and (|76p are the closed system for the numerical investigation 
of initial value problems such as the pair creation in strong field with the account of a 
back-reaction of the produced particles on the evolution of the mean gluon field. 

The simple solutions of the type (J7T|) do not satisfy to the Eq.((H2J) for the SU(3) case. 
But the reduction to the ordinary differential equation is possible, nevertheless, at use the 
more complicated representation of the type 



W a (t,p)=rj a W{t,p), a = 1,2 ... N. (77) 

We obtain the non-linear equations system for the admissible values of rf by substituting 
these relations in Eqs.()62|) 



a = 1,2... 8, (7£ 



where d abc has three independent non-zero values only 



7 Comparison with QED 

The equations (|35)) - (j42j) can be transformed to QED case by formally setting the color 
matrices t a equal to the unit one. Then the commutators of the gauge fields with the 
spinor components vanish, whereas the anti-commutators give a factor 2: 

d T a + gE a dl{p)a = 4^, 
d T ^+dUxK+gd a ± (p)(b ll E a -b^) = 0, 

d T c^ + d^(x)cr + gd1(p)(E a c^-F^c /3 /) = pta-mK, 

a T cf + d v ± (x)c>t-d f i(xy ± + 
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+gd^ip)(E a c^ + F^cl-F^c1) = ep^e^, 
d T S + dl(x)d ± a +gd < l(p)(SE a -F^d^) = 2me, 

d T d»-d»(x)d\\+gdl(p)(d»E a + F^d\\) = 2p^e a ^, 

d T e + gE a dl{p)e + 2mrf" = -2p^ s e af3& . (79) 

These equations are corresponds the formulae (4.9) - (4.16) of the work [3U] (part 2). The 
system (J79J) can be reduced for the simple field (|57|) to three scalar ordinary differential 
equations [2B]> which allows the simple numerical investigation. 



8 Summary 

We have derived the system of KE for description of quark- ant iquark plasma created 
from vacuum under action of a strong quasi-classical gluon field. The single-time Wigner 
function formalism allows in contrast to other approaches of this kind to formulate cor- 
rectly the Cauchy problem. It is particularly important for non-perturbative description 
of vacuum particle creation. 

We have analyzed some special cases of obtained system of KE (vacuum solution, 
space-homogeneous time dependent color electric field e.t.c). It is shown that KE is 
complex system of the partial differential equations even in the most simple case (chiral 
limit). That system is rather difficult for the numerical investigation, while the KE in 
QED allows the reduction to the set of ordinary differential equations [2E1 E2] for some 
field configurations. Thus the transition to QCD either sets higher request to level of 
computer calculations or needs a some additional non-perturbative model assumptions. 

Acknowledgments 

We thank A. Reichel for helpful discussions. This work was supported partly by the 
Ministry of Education of the Russian Federation under grant N E02-3. 3-210 and Russian 
Fund of Basic Research (RFBR) under grant 03-02-16877. 

References 

[1] U. Heinz and P.K. Kolb, |hep-ph/0204061| 



[2] S.K. Wong, Nuovo Chimento, A65, 689 (1970). 

[3] H.-Th. Elze and U. Heinz, Phys. Rep. 183, 81 (1989). 

[4] U. Heinz, Ann. Phys. 161, 48 (1985); 168, 148 (1986). 

[5] H.-Th. Elze, M. Gyulassy, and D. Vasak, Nucl. Phys. B 276, 706 (1986). 



14 



[6] M. Gyulassy and A. Iwasaki, Phys. Lett. B165, 157 (1985). 
[7] H.-Th. Elze, Z. Phys. C38, 211 (1988); C47, 647 (1990). 
[8] S. Mrowczynski, Phys. Rev. D39, 1940 (1989). 

[9] J.-P. Blaizot and E. Ianku, Nucl. Phys. B557, 183 (1999); Phys. Rep. 359, 355 
(2002). 

[10] G. C. Nayak, A. Dumitru, L. McLerran and W. Greiner, Nucl. Phys. A687, 457 
(2001). 

[11] A. K. Gangulu, P. K. Kaw, and J. C. Parikh , Phys.Rev. D48, R2983 (1993). 
[12] D.D. Dietrich, G. C. Nayak, and W. Greiner, Phys. Rev. D64, 074006 (2001). 
[13] G.C. Nayak and V. Ravishankar, Phys. Rev. D55, 6877 (1997); C58, 356 (1998). 
[14] A. Jain and V. Ravishankar, |hep-ph/021200l| 
[15] D.F. Litim and C. Manuel, Phys. Rept. 364 451 (2002). 

[16] Q. Wang, K. Redlich, H. Stocker and W. Greiner, Phys.Rev.Lett. 88, 132303 (2002); 
J.Phys. G28, 2115, (2002). 

[17] Yu.A. Markov and M.A. Markova, Transp. Theory Statist. Phys. 28 (1999) 645; J. 
Phys. G27, 1869 (2001); Phys.Rev. D64, 105009 (2001). 

[18] G. Gatoff, A.K. Kerman, and T. Matsui, Phys. Rev. 36, 114 (1987). 

[19] S.M. Schmidt, D. Blaschke, G. Ropke, S.A. Smolyansky, A.V. Prozorkevich, and 
V.D. Toneev, Int. J. Mod. Phys. E 7, 709 (1998). 

[20] S. Schmidt, D. Blaschke, G. Ropke, A. V. Prozorkevich, S. A. Smolyansky, and 
V. D. Toneev, Phys. Rev. D59, 094005 (1999). 

[21] J.C. Bloch, V.A. Mizerny, A.V. Prozorkevich, CD. Roberts, S.M. Schmidt, 
S.A. Smolyansky, and D.V. Vinnik, Phys. Rev. D 60 1160011 (1999). 

[22] S.A. Smolyansky, V.A. Mizerny, D.V. Vinnik, A.V. Prozorkevich, and V.D. Toneev, 
in "Progress in Nonequilibrium Green's functions", M. Bonitz (Ed.), World Scientific 
Publ., Singapore 2000; Proceedings of the Conference " Kadanoff-Baym Equations - 
Progress and Perspectives for Many-Body physics", p. 375, Rostock, 20-24.09.1999. 

[23] Y. Kluger, J.M. Eisenberg, and B. Svetitsky, Int. J. Mod. Phys. E2, 333 (1993). 

[24] Y. Kluger, E. Mottola, and J. Eisenberg, Phys. Rev. D58, 5015 (1998). 



15 



[25] D.V. Vinnik, R. Alkofer, S.M. Schmidt, S.A. Smolyansky, V.V. Skokov, and A.V. Pro- 
zorkevich, Few-Body System, 32, 23 (2002). 



[26 

[27; 

[28 
[29 
[30 

[31 

[32 

[33 
[34 



V.V. Skokov, S.A. Smolyansky, and V.D. Toneev, |hep-ph/0210"o"99"l 

D.V. Vinnik, A.V. Prozorkevich, S.A. Smolyansky, V.D. Toneev, M.B. Hecht, CD. 
Roberts, and S.M. Schmidt, Eur. Phys. J. C22, 341 (2001). 

I. Bialynicki-Birula, P. Gornicki, and J. Rafelski, Phys. Rev. D44, 1825 (1991). 

S. Ochs and U. Heinz, Ann. Phys. 266, 351 (1998). 

A. H611, V.G. Morozov, and G. Ropke, Theor. Math. Phys. 131, 812 (2002); 
132, 1029 (2002); |o5Ian"t-ph/0208083l 

S.A. Smolyansky, A.V. Prozorkevich, G. Maino, and S.G. Mashnik, Ann. Phys. 277, 
193 (1999). 

A. A. Grib, S.G. Mamaev, and V.M. Mostepanenko, Vacuum Quantum Effects in 
Strong External Fields, (Friedmann Laboratory Publishing, St. Peterburg, 1994). 

R. Hakim, Ann. Phys. 139, 230 (1982). 

F.J. Yndurain, Quantum Chromodynamics, (Springer- Verlag, N.Y., 1983). 



16 



